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We present a microscopic calculation of the 6 He /3-decay into the ground state of 6 Li. To this end, 
we use chiral perturbation theory at next-to-next-to-next-to-leading order to describe the nuclear 
weak-currents. The nuclear wave functions are derived from the J-matrix inverse scattering nucleon- 
nucleon potential (JISP), and the Schrodinger equation is solved using the hyperspherical-harmonics 
expansion. Our calculation brings the theoretical decay-rate within 3% of the measured one. This 
success is attributed to the use of chiral-perturbation-theory based mesonic currents, whose contri- 
bution is qualitatively different compared to standard nuclear physics approach, where the use of 
meson exchange currents worsens the comparison to experiment. The inherent inconsistency in the 
use of the JISP potential together with chiral-perturbation-theory based is argued not to affect this 
conclusion, though a more detailed investigation is called for. We conclude that any suppression of 
the axial constant in nuclear matter is included in this description of the weak interaction in the 
nucleus. 

PACS numbers: 23.40.-s,11.40.-q, 31.15.xj, 21.45.-v 



I. INTRODUCTION 

/3-decay is the every-day reflection of weak-interaction 
in nuclei. As such, it provides an experimental window to 
the properties of the weak interaction at nuclear density. 

In particular, theoretical studies of /3-decay rates of 
nuclei have argued for a suppression of the axial coupling 
constant gA, from its vacuum value, as extracted from 
the lifetime of the neutron g A = 1.2695 ± 0.0029 [l|, to 
unity, i.e. gA = 1 According to a recent study, this 

suppression occurs gradually, as the mass of the nucleus 
grows, and fully utilized for A w 40 0]. 

The ramifications of this suppression are numerous. 
For example, to the understanding of astrophysical phe- 
nomena, such as neutron-star cooling and core collapse 
supernovae, whose dynamics is controlled by the weak 
interactions. It is of no surprise that the source of 
this suppression has been the target of many theoretical 
works, which have associated it with a partial restora- 
tion of chiral symmetry in finite densities, deficiencies 
in the inclusion of correlations between nucleons, loop- 
corrections to the axial current originating in nuclconic 
excitations and mesonic currents, or combination of the 
three 0,10, MM- 

An important assumption is hidden in these suppres- 
sion mechanism: if a full calculation of the weak interac- 
tion inside nuclei was possible from first principles, then 
the calculated decay-rates should agree with the experi- 
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mental ones. That is, if one could describe correctly the 
correlations between nucleons, and the weak interaction 
of an external probe with a nucleus, then one should re- 
cover the physical value of the axial constant. In order 
to do that, pertinent is to solve the nuclear problem from 
first-principles. Due to the strong correlations involved 
in the problem, a calculation of nuclear wave functions 
from the nucleonic degrees of freedom is at reach only for 
very light nuclei. 

The lightest nucleus that undergoes a /3-decay is the 
triton. However, the theory cannot be checked in the 
triton since its half-life is used to remove some freedom 
in the weak interaction of a lepton with a nucleus, as will 
be explained cxplictly later. The lightest nucleus that can 
provide a test to the theory is thus 6 He. 6 He (J' r = + ) 
is an unstable nucleus, which undergoes a (3 decay with 
a half-life ti/ 2 = 806.7 ± 1.5 msec to the ground state of 
6 Li (,V = 1+) [13. 

However, a microscopic calculation of 6 He from its nu- 
clconic degrees of freedom, failed to reproduce the (3- 
decay rate. This study, accomplished by Schiavilla and 
Wiringa [ll|, has used the realistic Argonne ul8 (AV18) 
nucleon-nucleon potential, combined with the Urbana-IX 
(UIX) three-nucleon-force (3NF), to derive the nuclear 
wave functions, through the variational Monte-Carlo ap- 
proach. The model used for the nuclear weak axial cur- 
rent includes one- and two-body operators. The two- 
body currents are phenomcnological, with the strength 
of the leading two-body term - associated with A-isobar 
excitation of the nucleon - adjusted to reproduce the 
Gamow- Teller matrix element in tritium /3-decay. The 
calculated half-life of 6 He overpredicts the measured one 
by about 9%. An unexpected result of the calculation, 
was that two-body currents lead to a 1.7% increase in 
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the value of the Gamow- Teller matrix element of 6 He, 
thus worsening the comparison with experiment. The 
authors of this paper have presumed that the origin of 
this discrepancy is either in the approximate character 
of the VMC wave functions, or in the discrepancies of 
the nuclear model of the weak interaction. Pervin et al. 
[IH, have used the GFMC approach to evolve the VMC 
wave functions ansatz. They showed that this brings the 
single nucleon Gamow- Teller matrix element to about 0- 
3% deviation from the experimental value. However, the 
MEC are still expected to increase the deviation from the 
experiment to about 2-5% from the experimental value, 
leaving this problem intact. 

In the current paper, we argue that the origin of the 
discrepancy is indeed in the model of the weak interac- 
tion inside the nucleus. The foundation of such an argu- 
ment has to be in the underlying theory, i.e. quantum 
chromodynamics (QCD). Thus, we describe the weak 
currents within the nucleus, using an effective theory 
of QCD, namely chiral perturbation theory (^PT), ap- 
plicable at low energies, relevant to /3-decay processes 
use the triton /3-decay to calibrate 
the strength of the contact interaction part of the meson- 
exchange currents, thus the calculation is without any 
free-parameters. The six-body nuclear problem is solved 
in a fully ab-initio approach, expanded in hypcrspherical- 
harmonics function, from its nucleonic degrees of free- 
dom [It], [H, Gil ■ The nuclear wave functions are derived 
from J-matrix inverse scattering nucleon-nucleon poten- 
tial (JISP), describing two-nucleon scattering data and 
bound and resonant states of light nuclei to high accu- 
racy 0, [21], [22| . Using this approach not only brings 
the calculated /3-decay rate to within 3% of the mea- 
sured data, but also changes qualitatively the contribu- 
tion of the two-body meson exchange currents (MEC) 
compared to the work of Schiavilla and Wiringa. %PT 
based MEC are found to decrease the Gamow- Teller ma- 
trix element, compared with the increase found by Schi- 
avilla and Wiringa. We argue that this qualitatively dif- 
ferent behavior originates in the use of x?T based MEC, 
rather than the specific choice of the potential. 



II. THEORETICAL FORMALISM 

We start with a brief reminder of /3-decay process, and 
the formalism used in the calculation. The decay is a 
weak process, in which an unstable nucleus of charge Z 
emits an electron and anti-electron-neutrino, leaving a 
nucleus of charge Z + 1. The interaction is mediated 
through the exchange of heavy W + boson. As the mo- 
mentum transfer in the process is much smaller than the 
mass of the W + boson, the weak interaction Hamilto- 
nian is given by Hw = — G ^ d ^ / d 3 xj~(x) J +IJ/ (x), where 

G = 1.166371(6) x lfr n MeV" 2 is the Fermi coupling 
constant Q, V ud = 0.9738(4) is the CKM matrix ele- 
ment mixing u and al quarks involved in the process [l[, 



j~(x) is the lepton charge lowering current, and J +M is 
the nuclear charge raising current. The decay rate can 
be calculated using Fermi's Golden rule, and it is propor- 
tional to the squared matrix element of this weak Hamil- 
tonian (f\\Hw\\'i), where i (/) is the initial (final) state. 
The lepton current is well-approximated as a current of 
charged Dirac particles, thus results in kinematical fac- 
tors to the decay-rate. The weak nuclear current can be 

written as: J +fl = ^- ( J Vfl + J A ^ , where r + is a Pauli 

matrix. J Vfl (J A ^) has a polar- (axial-) vector symme- 
try. Here, we will discuss either a triton decay, or 6 He 
decay, hence the transitions are constrained by a selec- 
tion rule on the angular-momentum change in the transi- 
tion: A J = 0, 1. Thus, a multipole decomposition of the 
nuclear current is helpful. Due to the small momentum- 
transfer only the lowest multipoles contribute, i.e. the 
J = 1 electric multipole of axial- vector symmetry E A , 
and in the case of triton also the J = coulomb multi- 
pole of polar- vector symmetry Cq . We explicitly checked 
that indeed the contribution of neglected multipoles to 



poles to 



the decay-rate of 6 He can be bounded by 1% 
The leading order contribution to the E A and Cq op- 
erators are proportional to the Gamow- Teller and Fermi 
operators, respectively. Thus, it is customary, when dis- 
cussing the experimental rates, to talk about the empiri- 
cal Fermi and Gamow- Teller matrix elements, instead of 
E A and Cq , using the relations 



-in 



2J; + 1 



<c3T> 



and 



GT = 



6^ (gf) 

2Ji + l gA 



(1) 



(2) 



Here, (Cg) = (f\\C^\\i), and similarly for (E A ), J 4 is 
the total angular momentum of the initial nucleus, and 
g A = 1.2695 ± 0.0029 is the axial constant Q. 

As discussed by Simpson [24| . and later revisited by 
Schiavilla et al. [ll|, [25|, the "comparative" half- life is 
related to the "empirical" GT and F operators thorough 



(fT 1/2 ) t 



K/(G 2 \V ud \ 2 ) 
f+ %g A \GT\* 



(3) 



Here, K = 2^\n2/ml (such that K/(G 2 \V ud \ 2 ) = 
6146.6 ± 0.6 sec), and f A /fv = 1.00529 [H[ accounts 
for the small difference in the statistical rate function 
between vector and axial- vector transitions. Putting 
the measured 6 He comparative half-life (fT 1 / 2 )t = 
812.8±3.7sec [n|, one extracts |GT( 6 He)| cxpt = 2.161 ± 
0.005. For triton, (fT 1/2 )t = 1129.6 ± 3sec [H, thus 
|GT( 3 H)| cxpt = 1.6560 ±0.0026 

In order to complete a calculation, we have to specify 
the detailed structure of the weak-current, and to calcu- 
late the nuclear wave functions. These will combine to 
produce the theoretical E A , which will be compared to 
the experimental ones above. 
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III. xPT WEAK CURRENTS IN THE NUCLEUS 

The main difference between the current work and pre- 
vious ones, is the physical origin of the currents. The last 
two decades of theoretical developments have provided 
us with an effective theory of QCD, in the form of xPT 
0, Q EE El- The xPT Lagrangian is constructed by 
integrating out QCD degrees of freedom of the order of 
A x ~ 1 GeV and higher. It retains all assumed symme- 
try principles, particularly the approximate chiral sym- 
metry of the underlying theory. This SU(2)a xSU(2)y 
symmetry is based on the small up- and down-quarks 
masses (compared to the QCD breaking scale). The lack 
of parity doublets in the QCD scale is interpreted as an 
indication that this symmetry is spontaneously broken, 
with the pions as the Goldstone-Nambu bosons. Their fi- 
nite, albeit small, mass is due to the finite quark masses, 
explicitly breaking the chiral symmetry. 

Furthermore, the chiral Lagrangian can be organized 
in terms of a perturbative expansion in positive powers of 
Q/A x where Q is the generic momentum in the nuclear 
process, i.e. the /3-decay or the pion mass [HI, [lj, [l5| . 
The Chiral symmetry dictates the operator structure of 
each term of the effective Lagrangian, however not the 
coupling constants. A theoretical evaluation of these co- 
efficients, or low-energy constants (LECs), is equivalent 
to solving QCD at low-energy, and it is not yet feasible to 
obtain them from lattice calculations because of compu- 
tational limitations. Alternatively, these undetermined 
constants can be constrained by low-energy experiments. 

As the chiral symmetry is a gauging of the electro- weak 
interaction, the weak currents are the Nother currents of 
this symmetry. The weak axial current adopted in this 
work is the Nother current derived from the axial sym- 
metry of the chiral Lagrangian up to N 3 LO [13, [H[. At 
leading order (LO) this current consists of the standard 
singlc-nuclcon part, which, as mentioned above, at low 
momentum transfer is proportional to the Gamow- Teller 
(GT) operator, 



E 



1 |LO : 



= «3a(67t) 



- 1/2 X>^ 

i=l 



(4) 



where ai, t- are spin and isospin-raising operators of the 
ith nucleon. 

Corrections to the single-nucleon current appear at 
N 2 LO in the form of relativistic terms. It is easily verified 
p7| that the single nucleon current achieved in the xPT 
formalism, is identical to that achieved in the standard 
nuclear physics approach (SNPA). 

At N 3 LO, additional corrections appear in the form of 
axial MEC. While the relativistic corrections are negligi- 
ble for the half life, the MEC have a substantial influence 
on this /3-decay rate. This is a reflection of the fact that 
Ei is a chirally unprotected operator (29|. The MEC, 
to this order, include two topologies: a one-(charged)- 
pion exchange, and a contact term (that represents, for 
example, two-pion exchange or the exchange of heavier 



mesons). In configuration space the one-pion exchange 
part of the axial MEC is given by: 

-?4 8 (ni) = OW A (r i i)+ 
9a 

+c 3 (T^' a - Te Q )TO 2 yJ A (ry )+ 
+|(0^-0*>KtfSA(r«)- 
-(c 4 + ^)ml (r® a yl K {rij) + |o® a 2/oA( r «) 

Where « 92.4 MeV is the pion-decay constant, M sa 
938.9 MeV is the mass of the nucleon, m, « 139.57 MeV 
is the charged-pion mass [l[, and the low-energy con- 
stants c 3 = —3.66(8) and 64 = 2.11(9) are calibrated in 
the 7r-N sector [3(| ■ The operators used here are defined 
as, 



0% 



__jL(f(i) x r^nP^ ■ fia + P 2 <? (1) • r ia ) 

2 Q — 1 ' 12' 12 j ) 

and © = x, + , — . In addition, the Yukawa- like functions 
are: 



2/oaW 



fi 3 h -, 1 



ViA(r) = 

7r / \ 1 did, 

2/2A r = — 2 r "5 — -E-VOA r • 
m% or r or 

S\ is a cutoff function, which we take as a Gaussian. 

Apart from this, the MEC include a contact term, that 
has the form: 

IMfi-jgim) = lo^sfini), (5) 

9A 



where the "smeared" delta function is 



(2tt> 



(6) 



The LEC d r is the only LEC up to N 3 LO that cannot 
be calibrated in the single nucleon sector, as it originates 
in the contact interaction 7T-NN in the chiral Lagrangian. 
As a result, in order to determine e? r , one has to resort to 
a larger nuclear system. We will use the triton half- life as 
an experimental datum to determine this LEC, Sec. IV Al 



IV. NUCLEAR WAVE FUNCTIONS 

The difference between the one-body contribution to 
the 6 Hc- 6 Li GT matrix clement and the experimental 
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value is of the order of few percent. A result which on 
the one hand is very satisfying, but on the other hand 
implies that numerical accuracy at a per mil level is re- 
quired if we to regard the 6 He /3-decay as a test of the 
MEC model. In view of this req uired level of convergence 
we use the JISP16 potential [21j| to model the interaction 
between the nucleons. The JISP16 NN potential utilizes 
the J-matrix inverse scattering technique to construct a 
soft nuclear potential, formulated in the harmonic oscilla- 
tor basis, that by construction reproduces the NN phase 
shifts up to pion threshold and the binding energies of 
the light nuclei with A < 4. 

We use the Hypersphcrical-Harmonics (HH) expansion 
to solve the Schrodingcr equation. The HH functions con- 
stitute a general basis for expanding the wave functions 
of an ^4-body system 3l|. In the HH method, the trans- 



TABLE I: The JISP16 NN interaction d He, d H binding en- 
ergies, rms matter radius, and the leading order GT matrix 
element as a function of K max . 



lational invariant wave-function is written as 









^He 






K 


B.E. 


radius 


B.E. 


radius 


GT lo 


4 


8.094 


1.632 


7.364 


1.653 


1.6656 


6 


8.233 


1.656 


7.512 


1.680 


1.6620 


8 


8.319 


1.677 


7.604 


1.704 


1.6575 


10 


8.351 


1.691 


7.641 


1.720 


1.6547 


12 


8.360 


1.697 


7.651 


1.727 


1.6538 


14 


8.365 


1.701 


7.657 


1.733 


1.6530 


16 


8.367 


1.704 


7.660 


1.736 


1.6526 


18 


8.367 


1.705 


7.661 


1.738 


1.6524 


[2Uv 


8.354 




7.648 






[21]e 


8.496(20) 




7.797(17) 






Exp. 


8.482 


7.718 





* = C n[K]Rn(p)y[K}(^,Si,ti) (7) 
n[K] 

where p is the hyperradius, and R n (p) are a complete set 
of basis functions. The hyperangle, ft, is a set of 'iA — 4 
angles, and IVW](fti s i, ti) are a complete set of antisym- 
metric basis functions in the Hilbert space of spin, isospin 
and hypcrangles. The hyperradius p is symmetric un- 
der particle permutations since p 2 = j( r i ~ r j) 2 - 
The functions y[K] (ft, Sj, U) are characterized by a set 
of quantum numbers [K] [l8l . [l9| and possess definite 
angular momentum, isospin, and parity quantum num- 
bers. They are the eigenfunctions of the hypersphcr- 
ical, or generalized angular momentum operator K 2 , 
K 2 y [K] (ft, Si , U) = K(K + 3 A- 5)y [K] (ft, Si , U). The de- 
tails of our method are explained thoroughly in Ref. [32| . 



V. RESULTS 
A. The triton /3-decay - Calibration of d r 

Our results for the ground state properties of the A = 3 
nuclei, 3 H and 3 He , are presented in tablefl] In the table, 
we present the energies, matter radii, and the leading 
order GT matrix element (see Eq. (H}) as a function of 
K ma x, the limiting value of the hypcrspherical angular 
momentum K in the HH expansion. As we are using the 
bare interaction our results are variational. 

From the table it is evident that an excellent conver- 
gence is achieved for the A = 3 nuclei. Our results indi- 
cate that the JISP16 potential leads to an undcrbinding 
of about 80kcV for the 3 Hc and 120kcV for the triton. 
Comparing our results with the NCSM results of Shi- 
rokov et al. [2 ill . w e see a nice agreement with their vari- 
ational results [2l|v but a discrepenc y o f about 130keV 
with their effective interaction results [2l|]£. It should be 
noted that the GT matrix element converges much faster 
then the matter radius. This property can be probably 
attributed to the fact that the GT is a medium-range 



operator, which is influenced by the asymptotic behav- 
ior of the wave function, described correctly using the 
hyperspherical functions. Comparing the JISP16 leading 
order GT matrix element with those of other potential 
models, see table [TT1 we observe that the JISP16 poten- 
tial model leads to an enhancement of the 1-body matrix 
element and it almost coincides with the experimental 
value. This property is found also for the UCOM po- 
tential, and might be a result of the minimization of the 
contribution of 3NF to the binding energy, which is in the 
essence of both these potentials. In general, one observes 
from the table that non-local potentials, such as Bonn or 
the N 3 LO potentials, tend to predict a value for the GT 
matrix clement which is closer to experiment than the 
local potentials. 



TABLE II: The dependence of the triton /3-decay leading or- 
der GT matrix-element on the potential model. 



Potential model 



AV184-3NF [33J 
Bonn+3NF [34] 
Nijm+3NF [35] 
N 3 LO+3NF [39] 
UCOM [40J 
JISP16 [This work! 



Expt. 



GT|lo 
1.598(2) 
1.621(2) 
1.605(2) 
1.622(2) 

1.65(1) 
1.6524(2) 
1.656(3) 



As explained in Sec. IIII1 we use the triton half-life as 
an experimental input to determine the LEC d r . That 
is, we use the trinuclei wave functions to evaluate the 
matrix-clement | ( 3 Hc| \E^\ | 3 H) |, of the Ef operator built 
from the x?T based weak-current, as a function of d r , 
for various high-energy-cutoff values. Using the experi- 
mentally derived value for this matrix element we get the 
following calibration for <i r (A x ): 



d r (A x = 500 Me V) 
d r (A x = 600 MeV) 
d r (A x = 800 MeV) 



0.583(27) t (38) flA 
0.625(25) t (35) SA 
0.673(23) t (33) 9A . 



(8) 
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The numbers in parenthesis denote uncertainties in the 
last digits. The first error is due to the uncertainty in 
the triton half-life, whereas the second one is due to un- 
certainty in gA (the numerical error is negligible). 

B. The 6 He- 6 Li Gamow- Teller matrix-element 

Turning now to the A = 6 case, we present in table Hill 
our results for the ground state properties of the 6 He, and 
Li nuclei. As evident in the table, at the value K max = 
14, which corresponds to about 2 — 3 • 10 6 basis states, the 
binding energies of the 6-body nuclei are obtained with 
an accuracy of few hundreds keV. 



TABLE III: The JISP16 NN interaction 6 He, 6 Li binding en- 
ergies, rms matter radii, and the leading order GT matrix 
element as a function of K max - 





"He 




b Li 






-Kmax 


B.E. 


radius 


B.E. 


radius 


GT lo 


4 


18.367 


1.840 


19.392 


1.859 


2.263 


6 


24.103 


1.902 


26.124 


1.909 


2.247 


8 


26.392 


1.979 


28.854 


1.984 


2.234 


10 


27.560 


2.051 


30.156 


2.051 


2.232 


12 


28.112 


2.112 


30.797 


2.110 


2.229 


14 


28.424 


2.165 


31.132 


2.160 


2.227 


oo 


28.70(13) 




31.46(5) 




2.225(2) 


|2U 


28.32(28) 


31.00(31) 




Exp. 


29.269 


2.18 


31.995 


2.09 


2.170 



Taking a closer look at the table, we find that the 
binding energies exhibit an exponential convergence. De- 
ploying this observation we extrapolate our results to 
the K rnax — y oo limit, using the formula E(K max ) = 
Eoo + Ae~ maa . Fitting the parameters E^, A, a to 
the entries of table IHII in the range K max > 6, 8, 10 we 
find a rather stable value for Eoo with variance of about 
50keV for 6 Li and 130keV for 6 He. The resulting bind- 
ing energies are 28.70MeV for 6 He and 31.46MeV for 6 Li. 
While these results are roughly 550kcV below the exper- 
imental values, the difference AE = 2.76McV between 
the binding energies of the two nuclei differs by merely 
34kcV from the experimental value, AE = 2.726MeV. 
In the last column of table IIIII we present our 6 He- 6 Li 
leading order GT transition matrix element, i.e. at the 
Tbody level. It can be seen that the convergence pattern 
of the matrix element is not regular. Extrapolating its 
value using the expression GT(K max ) = GT 00 +Be _ ' 3Kma * 
for K max > 0, we get GT^ = 2.225(2). The fits of the 
extrapolation formulae to the calculated values are pre- 
sented in Fig. [T] for the binding energies and in Fig. [2] for 
the GT matrix-element. 

The value GT = 2.225(2) we obtained for the JISP16 
potential, is in accordance with the values GT = 2.28 
for AV8'/TM'(99) and GT = 2.30 for AV8' obtained by 
Navratil and Ormand [36], GT = 2.28 for the N 3 LO 
NN-force of Navratil and Caurier [37[, GT = 2.25 for 
AV18/UIX of Schiavilla and Wiringa [ll[, and GT = 



2.16 - 2.21 for AV18/IL2 by Pervin et al. Q3. Moreover, 
it can be seen that our accuracy in estimating the GT 
matrix element is at the level of per mil. Such an accu- 
racy enables us to disentangle numerics from physics and 
validates the use of the 6 He /3-decay as a testing ground 
for an axial MEC model. 

Incorporating the %PT based contributions to the 
weak-current we can finally calculate the full 6 He- 6 Li 
GT matrix-element at the N 3 LO level. In table IIV1 we 
present the transition matrix-elements as a function of 
Kmax and the cutoff A x . The appropriate values of d r 
are taken from Eq. ([5]). Two important observations 
can be drawn from the table, (i) the numerical accuracy 
of the calculated GT matrix-element is few per mil, and 
(ii) there is only a very weak dependence on the cutoff 
A x , which is of the same order of magnitude. The sec- 
ond observation implies that there is no need to refine 
our calculation, and moreover, the contribution of higher 
order x?T corrections to the weak-current are negligible. 

TABLE IV: The dependence of the full (l-body+2-body) 6 He- 
6 Li GT matrix-element on K max as a function of the cutoff 
A x , at the N 3 LO level. 

K ma x A x = 500MeV A x = 600MeV A x = 800MeV ' 



4 


2.1870 


2.1798 


2.1703 


6 


2.1850 


2.1805 


2.1746 


8 


2.1868 


2.1850 


2.1826 


10 


2.1937 


2.1932 


2.1927 


12 


2.1951 


2.1952 


2.1955 


11 


2.1970 


2.1975 


2.1983 



Summarizing, the predicted GT of 6 He is: 

|GT( 6 He)| thco = 2.198(l) A (2) N (4) t (5) gA = 2.198 ± 0.007 

(9) 

The first error is the cutoff variation dependence, the 
second is numerical, the third is due to uncertainties in 
the triton half-life, and the last is due to uncertainties in 
gA- This should be compared to the experimental matrix- 
element |GT( 6 He)| cxp t = 2.161±0.005. Thus, the theory 
overpredicts GT by about 1.7%. 

VI. DISCUSSION 

The use of phenomenologically based potential, JISP, 
combined with a %PT based MEC, is an inconsistency in- 
herent to our calculation. Clearly, the chiral Lagrangian 
can be used to derive the nuclear forces as well. This 
inconsistency, however, allows us to overcome limited 
computational resources, as well as theoretical difficul- 
ties (the N 3 LO nuclear potential has not been fully de- 
veloped yet), and to accomplish the task of a microscopic 
calculation of a six-body problem. The use of a hybrid 
approach, sometimes coined EFT*, has had great suc- 
cess in the literature [27l [28l. |38|. In all these checks, the 
phcnomcnological nuclear forces included realistic poten- 
tials, the AV18+UIX force model. This potential, though 
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FIG. 2: (Color online) The convergence of the GT matrix element for the 6 He- 6 Li /3-decay. The continuous line is the fit 
GT(K max ) = GToo + Be _/3Kmax , the dashed line is the extrapolated value GT^. 



different than the X-PT force models in the short-range 
character of the force, has the correct long-range behav- 
ior, due to the pion exchange. The JISP potential, how- 
ever, is different in this respect, as it is built in an ab-exitu 
approach, and does not have an asymptotic long pion be- 
havior, thus not consistent with chiral symmetry even at 
long-distances. In addition, the JISP potential does not 
include a three-body force. 

It is hard to estimate the effect of these approxima- 
tions. However, in a recent work [39], the triton /3-decay 
process was calculated using force model and current de- 
rived consistently from the same xPT N 3 LO Lagrangian. 
One of the conclusions of this work has been that the 
short-range correlations of the force and the short-range 



correlations of the weak current are not correlated, thus 
the effect of the three-body-force is negligible for GT-type 
operators. In addition, the JISP potential successfully re- 
produces nucleon-nucleon scattering data, and the bind- 
ing energies of A < 16 mass nuclei. However, the most 
convincing reason to believe the stability of the current 
results, is the minimal dependence of the half-life in the 
cutoff. 

We thus believe that even in the current calculation, 
the effect of the approximation will not change qualita- 
tively the results, and the effect of the MEC. The qual- 
itative difference originates in the different structure of 
the SNPA and X PT based MEC. 

A careful analysis of the difference between the MEC 
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originating in xPT and those used in SNPA, has been 
accomplished by Park et al [27| . They have shown that 
one-pion exchange term exists in both models. Of par- 
ticular importance is the part of this term in the SNPA 
based MEC that represents the exchange of a pion due to 
a delta excitation of the nucleon, which is found to cor- 
respond roughly to the £3 term in the xPT based MEC. 
The coupling constant of this term g^NA has been fixed 
by Schiavilla and Wiringa [h[ , so that the theory would 
reproduce the triton half-life. 

However, differences between the approaches arise in 
their short range character. In the SNPA approach, these 
correspond to the exchange of a p-meson. Such a term 
does not exist in the xPT approach as it arises only at 
N 5 LO " 



27] 



Moreover, a contact interaction of the form 
of Eq. (J6]) does not appear in the SNPA approach. It is 
this contact interaction that creates the qualitative dif- 
ference between the current work and that of Schiavilla 
and Wiringa [ill ]. 
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FIG. 3: (Color online) Relative contributions to the theoret- 
ical GT matrix elements as a function of the EFT cutoff. All 
the results are normalized to the empirical values. The blue 
lines indicated by A = 3 correspond to the 3 H- 3 He /3-decay. 
The red lines indicated by A — 6 correspond to the 6 He- 6 Li 
case. Dashed lines correspond to the 1-body impulse approx- 
imation (IB). Dashed-dotted lines correspond to 1-body plus 
one-pion exchange current (OPEC). Continuous lines corre- 
spond to full calculation (note that in the case of 3 H this is 
calibrated to give exactly the experimental value). 



In order to acknowledge that, we plot in Fig. [3] the 
relative contribution of each the terms, i.e. one-body, 
one-pion-exchange and full calculation, to the GT matrix 
element. One first recognizes that the one-pion contribu- 
tion to the matrix element has a positive sign in both 3 H 
and 6 He, and that the contact interaction has a negative 
contribution to the matrix element. In the case of 3 H 
this is only a partial cancellation, as it is calibrated to 
increase the 1-body matrix element and to bring the cal- 



culation into the experimental value. In view of the fact 
that the one-body calculation in the case of the JISP po- 
tential almost exhausts the total GT strength, one might 
suggest that the negative sign of the contact interaction, 
as well as the partial cancellation is an artifact of the po- 
tential. However, the same partial cancellation is found 
also in a consistent N 3 LO calculation of 3 H decay, thus 
it is not a result of the use of the JISP potential [39| • 

In contrast to 3 H, when examining the case of 6 Hc, 
one observes that the negative contribution of the con- 
tact term is bigger (in absolute value) than the one-pion- 
exchange contribution, thus leading to a total negative 
contribution of the MEC. This negative contribution is 
needed as the single-nucleon GT is bigger than the ex- 
perimental GT. 

Recalling the fact that the SNPA approach does not 
contain a contact interaction, we understand the origin 
of the positive contribution of the MEC in that approach, 
which increases the difference between the calculated and 
measured decay-rates. 



VII. SUMMARY AND CONCLUSIONS 

In this work we have used the 6 He /3-decay as a testing 
ground for the nuclear weak-current derived from xPT. 
A precondition for such a task is an accurate evaluation 
of the 6 He- 6 Li weak transition matrix-element at the per 
mil level. To this end we have used the soft NN potential 
JISP 16 to describe the nuclear dynamics and the HH ex- 
pansion method to solve the Schrodinger equation. The 
weak interaction in the nucleus is completely determined 
by fixing the short range behavior of the scattering oper- 
ator to reproduce the experimental 3 H half-life, resulting 
in a parameter-free prediction of the 6 Hc /3-decay rate. 

We have found that at the 1-body, impulse approxi- 
mation, level the 6 He- 6 Li GT matrix-element is over pre- 
dicted by roughly 3%. This observation for the JISP 16 
potential is in agreement with previous findings for other 
potential models. Adding 2-body, meson-exchange, cur- 
rents derived within xPT, we have found that in contrast 
with the previous work of Schiavilla and Wiringa [ 1 II ] , the 
2-body MEC contribution to the 6 He- 6 Li transition ma- 
trix element is negative. We argue that this difference 
originates in the different short-range character of the 
MEC derived in the two approaches. We find that both 
for 3 H and 6 He, there is a sign difference between the 
positive contribution of the long-range one-pion-exchange 
current, and the negative contribution of the contact in- 
teraction in xPT, representing higher degrees of freedom 
which were integrated out in the development of the effec- 
tive theory. In the case of the 6-body transition, however, 
the contact interaction has a bigger value than the one- 
pion exchange contribution. Thus, it provides the origin 
to the sign difference between the MEC contribution in 
3 H and 6 He. This contact interaction does not exist in 
the standard nuclear physics approach, adopted by Schi- 
avilla and Wiringa. Therefore, the reconciliation between 
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the theoretical and the experimental 6 He half- life is due 
to the use of the x?T formalism. 

Our calculation points to an agreement at the level 
of about 1.7% between the measured and calculated GT 
matrix elements. This result should be contrasted with 
the difference of 5.4% obtained by Schiavilla et al. [Til ] . 
More importantly, it shows that dominant contributions 
that arise naturally in the xPT formalism, and do not 
appear in the standard nuclear physics approach, are es- 
sential to a successful prediction of this weak observable. 
In order to pin-point this argument, a use of a consistent 
approach, in which both the weak currents and the nu- 
clear forces are derived from the same microscopic theory, 
is called for. 

The agreement between the calculated and measured 



decay rates of 6 Hc indicates that there is no signature 
in this observable for an additional suppression of the 
axial constant. It appears that all the needed suppres- 
sion originates in correlations between nucleons in the 
nucleus, revealing itself in the form of exchange currents. 
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